Non-universality of the scaling exponents 
of a passive scalar convected by a random flow 
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We consider passive scalar convected by multi-scale random velocity field with short yet finite 
temporal correlations. Taking Kraichnan's limit of a white Gaussian velocity as a zero approximation 
we develop perturbation theory with respect to a small correlation time and small non-Gaussianity 
of the velocity. We derive the renormalization (due to temporal correlations and non-Gaussianity) 
of the operator of turbulent diffusion. That allows us to calculate the respective corrections to the 
anomalous scaling exponents of the scalar field and show that they continuously depend on velocity 
correlation time and the degree of non-Gaussianity. The scalar exponents are thus non universal as 
was predicted by Shraiman and Siggia on a phenomenological ground (CRAS 321, 279, 1995). 
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The most striking feature of turbulence is its intermit- 
tent spatial and temporal behavior. Statistically, inter- 
mittency means substantial non-Gaussianity. For devel- 
oped turbulence, where the correlation functions are scale 
invariant at the inertial interval of scales, the intermit- 
tency is manifested as an anomalous scaling of correlation 
functions. That means that some random field 8(r, t) has 
the structure functions S2n — ([&{t, r l) — 0(t, T2)] 2n ) oc r^™ 
with the exponents (2-n that are not equal to nC,2- As a 
result, the degree of non-Gaussianity, that may be char- 
acterized by the ratio ^n/S'J , depends on scale. Exper- 
iments and simulations show that the intermittency and 
anomalous scaling of the scalar field passively convected 
by a fluid are much stronger pronounced than the in- 
termittency of the velocity field itself It is in the 
problem of passive scalar where consistent analytic the- 
ory of an anomalous scaling starts to appear ^ 
is intuitively clear that the physical reason for scalar in- 
termittency is a spatial inhomogeneity of the advecting 
velocity. The analysis of the velocity field with smooth in- 
homogeneity shows, however, that there is no anomalous 
scaling of the scalar whatever be the (finite) temporal 
correlations of the velocity ||||,|l3|] . Analytic treatment 
of a non-smooth multi-scale velocity was possible hith- 
erto only in the so-called Kraichnan's problem of white 
advected scalar jj| where the correlation functions sat- 
isfy closed linear equations of the second order Q . It has 
been shown that, even without any temporal corre- 
lations, spatial non-smoothness of the velocity provides 
for an anomalous scaling of the scalar. The anomalous 
parts appeared as zero modes of the operator of turbu- 
lent diffusion and entered the correlation functions due 
to matching conditions at the pumping scale The 
coefficients at the modes were thus pumping-dependent 
while the form of any zero mode was universal i.e. deter- 
mined only by the exponent of the velocity spectrum and 
space dimensionality. In particular, the exponents £ n of 
the scalar were universal for the delta-correlated velocity. 



Now, what is the role of velocity temporal behavior in 
building up intermittency of the scalar field? It was ar- 
gued phenomenologically by Shraiman and Siggia [0 that 
the exponents of the scalar field depend on more details 
of the velocity statistics "than just exponents" . Here, we 
consider the simplest possible generalization of Kraich- 
nan's problem and consistently derive the equations for 
scalar correlation functions in the case of short yet finite 
velocity correlation time r r which is supposed to be a 
power function of the scale r. The behavior of the ratio 
r r /t r is important, where t r is the turnover time at the 
scale r. If the ratio tends to zero at decreasing r then we 
return to the (^-correlated case. If the ratio increases at 
decreasing r we encounter the problem of the quenched 
disorder type which should be considered separately. We 
consider the marginal case of a complete self-similarity 
where e = T r /t r does not depend of r and formulate the 
perturbation theory regarding the ratio as the small pa- 
rameter of our theory. We show that £2 does not depend 
on e while £ n for n > 2 are e-dependent that is the set of 
the exponents is non universal along with the prediction 
of Q. The principal difference between the second and 
higher correlation functions is naturally explained on the 
language of zero modes: there is no zero mode (except 
constant) for the pair correlator while the zero modes 
of the high correlators depend on the precise form of 
the operator of turbulent diffusion which is e-dependent. 
This is formally similar to what has been discovered by 
Kadanoff, Wegner and Polyakov in the theory of phase 
transitions: the critical exponents continuously depend 
on the amplitude of the operator term with dimension d 
added to the Hamiltonian 

Note that the below results cannot be directly applied 
to the description of scalar advection by a Kolmogorov 
turbulence: because of sweeping effect, the different-time 
velocity statistics is not scale-invariant in the Eulerian 
frame AO . Our use of scale- invariant velocity is intended 
to establish the general fact of the sensitive dependence 
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of scalar exponents on the velocity statistics. 

The advection of passive scalar 6(t, r) by an incom- 
pressible flow is governed by the equations 

{d t -P)9 = <j>, P(t) = ^u a V a + KA, V a v a =0, (1) 

where n is the coefficient of molecular diffusion. The 
advecting velocity v and the source are independent 
random functions. A formal solution of (Q) is 

0(t,r)=J (ftiTexpQ* di'P(i')^(*i,r), (2) 

where Texp designates the chronologically ordered ex- 
ponent. From (|^) it follows 

F„(t,r 1 ,...,r 2 „) = (8{t,r 1 )...d(t,r 2n )) 

t t . „„ . 



A)o(t) = «WV>f(0)Vf). 



(8) 



= Jdh... Jdt 2n A/f[<t>{U,ri)\ , 

— oo — oo \i=t I 

A=(Q), Q(t) = JJ Texp (j( <P«,r 2 



(3) 
(4) 



Differentiating .4 over the current time i, one gets 

M = (v(t)Q(t)\ , P = - «?V? , (5) 



where Vj = v(i, rj) and V» = djdvi. Here and below 
summation over both repeated vector superscripts and 
subscripts enumerating points is implied. The identity 
(||) can be brought to the form 



d t A(t) = KVlA{t) + / dt'Af(t')A{t - t') , 



(0) 



where TV (i) is to be found. The decay of Af(t) is deter- 
mined by the velocity correlation time r r which is sup- 
posed to be much smaller than the spectral transfer time 
characteristic of (^). It is the reason why the upper limit 
in (JsJ) can be substituted by infinity. 

We shall find the first terms of the expansion of TV in 
T r . Let us first examine the Gaussian contribution to dtA 
related to reducible correlation functions of v 



dtvf (t) V? ( T exp / dt'V{t') v$ (t)Vj(5(t) 



(7) 



where the product vf(t)Vj(i) should be substituted by 
the pair correlation function (v a (t, ri)vP(t, rj)). The 
integrand is nonzero for t — t < r r and consequently 



Texp 



fi dt'V(t') can be expanded over t — t. The zero 



The first and the second terms of the expansion produce 
the linear in r r contribution 



t t 
TVi(t) = J dt x 
o 



*afe?(*)V?«2(ti)V24(t 2 )Vj!%(0)Vg, 
«?(*)V?!>2(ti)V2C(«a)V^(0)Vf] (9) 



-«ti^(ti)v?vK(o)vf 



where the products vv and vv should be substituted by 
the corresponding pair correlation functions. 

For a short-correlated velocity field, the leading non- 
Gaussian contribution to the correlation functions of v 
is determined by the irreducible part of the fourth-order 
correlation function of v. Generalizing the trick leading 
from (||) to (R) we obtain the non-Gaussian term TV" n G(i): 



dtx rf«2((^(t)V?<(t 1 )V^(t 2 )Vf<(0)V^)>, (10) 



where double angular brackets stand for the cumulant. 
The operator A(t) is exponential in time 



A{t) = exp[(t - t ){ K V* + £)]A(t ) 



(11) 



asymptotically at t — t ^> r T . Substituting ( |ll| ) into (|£|), 
expanding expit'C) and keeping only the principal terms 
we find the operator of turbulent diffusion 



C — Co + C\ + £j + C n c 



term gives 



(12) 



£{o,i,„G} = / dtU {0XnG} (t), £[=- dttN {t)L . 
Jo Jo 

Using (||,[ll]) we can obtain the expression for dtF n . For 
the pumping (5-correlated in time, one gets 

d t F n (t, n, . . . ,r 2n ) - CF n (t,Vi, . . . ,r 2 „) 

= M[xi2F n -i(t, r 3 , . . . , r 2 „) + permutations] . (13) 

Here, the function x( r i2) = / dt {<j>(t, ri)^>(0, r 2 )} decays 
on the pumping scale L and x(0) is the production rate of 
6 2 . The operator M. in (|l3| ) can be estimated as A(tl)- 
The account of temporal correlations of the pumping 
(which can be done perturbatively as long as the pump- 
ing correlation time is much less than the time of scalar 
transfer) results in an extra renormalization of M. oper- 
ator. Its explicit form is unimportant for what follows. 
Indeed, the balance equation ( |l3| ) contains the renormal- 
ization (due to velocity temporal correlations and non- 
Gaussianity) of all three relevant quantities: pumping, 
turbulent diffusion and molecular diffusion (the last term 
in TVi). We discuss here only the scaling exponents in the 
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convective interval of scales (see below) that are deter- 
mined solely by the form of the operator of turbulent 
diffusion C. 

Let us consider the pair correlation function of the ve- 
locity to be scale-invariant: 



([v a (t, r)-v a (0, 0)K(i, r)-«"(0, 0)}) = 2K^(t, r), 

£> r 2-7 



K 



a/3 _ ; 



•a/3 



(14) 



with the correlation time r r — TL(r/L) z . Dimension- 
less functions g± and gu satisfy the incompressibility 
condition [d — i)g±(x) — zx a d[x 1 ~ a g\\ (x)]/dx where a — 
(2 — ~f)/z — 1. Their normalization is fixed by below ex- 
pressions (00,03). The main term in ((L2|) is m 



oo 

C = -J^lCf (ry)V?Vf Kf = 2f dtlC°0(t) (15) 
»j o 

a+4=5E *s£* or vr vf - \ j2 Bgvy] - 

/ dt 2 K a ?(t 2 ;r) / dt 3 K^(t 3 ;r) 
.Jti Jti 



]Cf = Dr 2 -< 



r a r l3 
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(16) 



The expressions ( |15| , |16|) lead to the following turnover 
time t r — (2 — -f)r n '/D'yd(d—l) obtained for the delta- 
correlated case [^jj?]]. Our marginal case corresponds to 
z = 7 and the small parameter of the perturbation theory 
is thus 



e = Dr L L-~< 1 d{d - l)/(2 - 7 ) < 1. 



(17) 



Note that e contains d 2 which tells us that the space di- 
mensionality should not be very large for the approxima- 
tion of a short correlation to be valid: the characteristic 
time of the scalar transfer (proportional to d~ 2 ) should 
be larger than the correlation time. 

Starting from the expression for the pair velocity cor- 
relator ( fl4| ) we can obtain the first Gaussian e-correction 
to (|l5|). Calculating ([)]) and then integrals in (12) we 
find 



K \ - 



pOO 

-2 / dtxV?VP 



2 ^ 



• OO POO 

dt 2 K^(t 2 ;r) / dhK 0v {h;v) 

1 Jt x 



(18) 



where = V?K%P, K^ v = V^K^ and 



OO 



B af3 (r) = eDr 
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b\i5 af3 + b ± ( S af} 



(19) 



Now we can analyze the equation (|T^) for F n . At the 
convective interval of scales L ^> r ^> [k(2 — ^) / D(d — 
l)] 1 /(2-^r) j the molecular diffusion term can be dropped: 
it is enough to require F n = at r = ||j7]]. Here, the 
zero modes of C are responsible for the anomalous scaling 
of F n ■ The scaling exponents of the bare operator Co and 
the perturbation operator C\ + C[ coincide. For a self- 
similar velocity statistics, the non-Gaussian contribution 
C n Q has the same scaling too. The first consequence is 

C ,o=dY,rlpd nj , C.-.j ^/-;, "(/•,,<•/;, -./),,) r 
e<f(2- 7 ) N - ;j , _ , .„,,„ 



that the exponent of the pair correlation function is 7 at 
arbitrary finite order in e for any 7 and d. Indeed, there 
is no zero mode of the two-point C with a nonzero posi- 
tive exponent that could provide an anomaly. Contrary, 
for n > 2, the account of the e-contributions to the bare 
operator Co should produce obviously e-dependent cor- 
rections to the exponents of zero modes and consequently 
e-dependent anomalous scaling. 

To illustrate the above conclusion about r-dependence 
of the scalar exponents, let us give an example where the 
calculation can be done explicitly. We consider a large di- 
mensionality (the limit 7 3> (2 — j)/d solved in 14 for 
r = 0) while assuming that, in addition to (0), l/d>e 
(it will be seen below how the parameter 7 enters the 
condition). The leading (in d) terms of the bare and 
the Gaussian perturbative operators in terms of relative 
distances 77, are as follows [multiplied by (2 — j)/dD] 



-8- B 



+ 46 ii 0) E 
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<1 - 2) 
87 



J2 r kidki+J-l+2b\ 

M>1 



i>3 



f '1:1 > f 



'TikU 



„2-7 r ip r jq 

3 r ip r jq 



d w d 3q + r\pr kl d t] du + (l6 - 4 7 - 8 (b™ +bf b^)) ]Tr^ 7 d 4: 



i>j,k>l 



i>j 



(20) 
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The summation is performed over n(n — l)/2 distances 
which are independent variables if d > n — 2. For the 
chosen form of K af3 {t,r), b\\ -> (2 - j){b\° ] d 3 + b^d 2 ), 

b± — ► (2 — 7)6^1°'' <i 2 at <i — > 00, with d-independent con- 



(0 



First, we calculate the corrections to the 



stants b j 

exponents related to the Gaussian correction (20) and 
then discuss the corrections due to non-Gaussianity. 

Solving the equation for the pair correlation function 
one can check that (2 = 7 is independent of e and d. 
Then, we consider the four-point correlation function 
(i, j, k, I = 1 . . .4). To get the main contribution at r <C L 
one has to perturb the bare zero mode of £0.0 + £0.1 ■ I n 
the limit under consideration, it is enough to consider 
only the mode 



= E(' 



"lif 1/2 £ (r. 



7 \2 
ik) 



(21) 



with the leading exponent A4(0) = 4(2— j)/d found in Jt]] 
by 1/d-expansion. The first e-correction to (pi] ) can be 
obtained by applying the operator — £q J(£i,ojCq o£o,i + 
£o,i£(7o^i,o) t° the mode Zo. The correction to the 
anomalous exponent is determined by the coefficient at 
ln(L/r) in the first e-contribution to Zo, it is 



A 4 (e) = A 4 (0) 



e(2-7) 
c?7 



4 + 67 - 2 7 2 



(22) 



The sign of the correction is positive for < 7 < 2. We 
can also calculate r-rclatcd corrections for the high-order 
functions by using the technique developed in juj for 
finding the largest anomalous exponents. For n-^^yd, 



A 2n (e) = n(n - l)A 4 (e)/2. 



(23) 



The scaling exponents thus depend not only on purely 
dimensionless quantities 7 and d yet also on a dimen- 
sionless ratio of dimensional quantities. In other words, 
the exponents depend on the form of the structural func- 
tions g± and gti. 

Considering opposite hierarchy of small parameters 
1/rf-Ce and neglecting 1/d corrections, one finds ^2n = 
71(2 at any order in e: temporal correlations by them- 
selves do not produce an anomalous scaling if it is absent 
in the uncorrelated case. In this limit, the anomalous 
exponents appear only in the next 1 jd order and are pro- 
portional to ed - see ju| for the details. 

Now let us discuss the non-Gaussian contribution to 
the anomalous exponents. We denote by e 4 the ratio 



of the cumulant to the fourth-order correlator and con- 
sider the limit 1 ~S> > e 4 d 3 . Using the expressions 
(|l^,^2|,^3|) we conclude that the contribution to A„ is 
proportional to n(n — l)e 4 c? 2 for n <C jd. 

To conclude, we learned that the scalar exponents are 
sensitive to the details of the velocity statistics. The ex- 
istence of two different contributions (due to temporal 
correlations and non-Gaussianity of the velocity) makes 
it possible that there exist some classes of the statistics 
with special relations between the contributions (some 
remarkable cancelations, for instance), their analysis is 
left for future studies. Hopefully, real turbulent flows 
belong to those classes and analytic expressions for the 
passive scalar exponents can be found some day. 
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